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THE UNIFIED METHOD FOR THE THREE- WAVE 
EQUATION ON THE HALF-LINE 
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Abstract. We present a Riemaim-Hilbert problem formalism for 

, . the initial-boundary value problem of the three-wave equation: 

Mh, bi — bj <r— » bk — bj bi — bk 
<f . Pij,t PH,x+yX WikPkj = 0, i,j,k = 1,2,3, 



on the half-line. 
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CO 
pH I 1. Introduction 

The 3-wave resonant interaction model described by the equations 

<N : Pij ,t - S-p^ + £ fc (^ - *e£W« = o, (i i} 

i,3,k=l, 2, 3; <n ^ aj, h ^ bj, for i ^ j, 

^2 ■ is one of the important nonlinear models with numerous applications to 

tj- ! physics [H] . The 3- and iV-wave interaction models describe a special 

^ , class of wave-wave interactions that are not sensitive on the physical 

nature of the waves and bear an universal character. This explains why 
they find numerous applications in physics and attract the attention of 
the scientific community over the last few decades [231 I2H [251 [201 HI 
1221 I2"T1 E] and the references therin. 

The 3-wave equations can be solved through the inverse scattering 
method due to the fact that equation ( II .ip admits a Lax representa- 
tion [211 [25]. But until the 1990s, the inverse scattering method was 
pursued almost entirely for pure initial value problems. In 1997, Fokas 
announced a new unified approach for the analysis of inital-boundary 
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2 J.XU AND E.FAN 

value problems for linear and nonlinear integrable PDEs [U EJ |3] . The 
Fokas method provides a generalization of the inverse scattering for- 
malism from initial value to IBV problems, and over the last fifteen 
years, this method has been used to analyze boundary value problems 
for several of the most important integrable equations with 2x2 Lax 
pairs, such as KdV, Schrodinger, sine-Gordon, and stationary axisym- 
metric Einstein equations, see e.g. [HE]. Just like the 1ST on the line, 
the unified method yields an expression for the solution of an initial- 
boundary value problem with that of a Riemann-Hilbert problem. In 
particular, the asymptotic behavior of the solution can be analyzed 
in an effective way by employing the Riemann-Hilbert problem and 
the steepest descent method introduced by Deift and Zhou [TO] . Re- 
cently, Lenells develop a methodology for analyzing initial-boundary 
value problems for integrable evolution equations with Lax pairs in- 
volving 3x3 matrices [7J. He also used this method to analyze the 
Degasperis-Procesi equation in [8]. 

Pelloni and Pinotsis also studied the boundary value problem of the 
N— wave equation by using the unified method [11]. Recently, Gerd- 
jikov and Grahovski considered Cauchy problem of the 3- wave equation 
with with non- vanishing initial values [12]. In this paper we analyze the 
initial-boundary value problem of the three-wave equation ( II. ip on the 
half-line. Compared with these two papers, there are two differences 
in our paper. The first difference is that we get the residue conditions 
of matrix function M in the Riemann-Hilbert problem (see (I2.27J) in 
the next section 2). The second difference is that we the jump matrix 
J is explicitly constructed ( see the equations ( I2.14p and ( I2.22p in the 
next section 2). Of course, the initial-boundary value problem for the 
3— wave equation does not need to analysis the global relation, because 
the initial data and the boundary data are all known. 

The organization of the paper is as follows. In the following section 
2, we perform the spectral analysis of the associated Lax pair. And we 
formulate the main Riemann-Hilbert problem in section 3. 
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2. Spectral Analysis 

Our goal in this section is to define analytic eigenfunctions of the Lax 
pair (12. ip which are suitable foe the formulation of a Riemann-Hilbert 
problem. 

2.1. Lax pair. We first consider the three-wave equations ( 11. II) . with 
(x,t) G Q, and Q denoting the half-line domain 

Q = {0 <x < +cx),0 < t <T} 

and T > being a fixed final time. We denote the initial and boundary 
values by Pij,o(x) and qij,o(t), respectively 

Pij,o 0) = Pij {%, 0) , q ijfi (t) = pij (0 , t) 

with Pijft(x) and qijfiif) are rapidly decaying. Equation (II. ip admits 
the following Lax representation [L 



(2.1) 



4>x = M<j>, 
<t>t = N(f>, 
where M = iXA + P and TV = iXB + Q, with 



A 



B 



( ai \ 

a 2 

V o 3 / 

/ &i \ 

fr 2 

V o o b 3 J 



(2.2) 



(0 Pl2 Pl3 ^ 
Vl\ p 2 3 
P31 P32 / 

/ n X2 P\2 n 13 p 13 \ 

Q = n 2 lP21 U23P23 

\ ^31^31 n 32 p 3 2 / 

Obviously, trace(A) = trace(B) = 0. We also assume that a\ > (12 > 
a 3 and b\ < b2 < b 3 . By using transformation 



// 



-iXAx-i\Bt 



we change Lax pair (2.1) in the form 



fi x - i\[A, fj] = P/j,, 
\i t - iX[B,fi} = Qfi, 



(2.3) 
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which can be further written in differential form as 
d(e~ i)Jx - ix6t fj,) = W(x,t,X), 

where 

W(x, t, A) = e~ ixAx - iX§t (Pdx + Qdt)fi 

and e A X = e A Xe~ A . 



(2.4) 
(2.5) 



2.2. Spectral functions. We define three eigenf unctions {^j}l of (12. 3p 
by the Volterra integral equations 

tij(x,t,\)=I+ [ e {ixAx+ix6t) Wj{x',i!,X). j = 1,2,3. (2.6) 



ij 



where Wj is given by (12. 5p with fj, replaced with \x^ and the contours 
{7j}? are showed in Figure 1. And we have the following inequalities 












.(x,i) 



.(x,i) 



(ar,t) 



O 



o 



o 



7i 



72 



73 



Figure 1. The three contours 71,72 and 73 in the (x, t)— domain. 



on the contours: 



7i 

72 
73 



x-x' > 0,t-f < 0, 
x-x' > 0,t-t' > 0, 
x — x' < 0. 



(2.7) 



So, these inequalities imply that the functions {^j}\ are bounded and 
analytic for A G C such that A belongs to 



/'1 

/'2 
/'.3 



(£2,0, a; 



(2.* 



ODi,0, £2), 

where {D n }\ denote two open, pairwisely disjoint subsets of the com- 
plex A plane showed in Figure 2. 
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D 2 

Figure 2. The sets D n , n = 1,2, which decompose the 
complex A;— plane. 

And the sets {D n }\ has the following properties: 

D 1 = {A G C|Re/i < Re/ 2 < Re/ 3 , Rez x > Rez 2 > Rez 3 }, 
D 2 = {X E C|Re/i > Re/ 2 > Re/ 3 , Rez 1 < Rez 2 < Rez 3 }, 

where k(X) and Zi(X) are the diagonal entries of matrices iXA and iXB, 
respectively. 

2.3. Matrix valued FUNCTIONS M n 's. For each n = 1,2, de- 
fine a solution M n (x,t, X) of (12. 3 j) by the following system of integral 
equations: 



(M n )o-(x,t, A) = <%+ / (e(- aAr - iABt ^ n (x',t', A)) y , A G D n , i, j = 1,2,3. 

(2-9) 
where PFn an M n are given by (12. 5p . and the contours 7™-, n = 1,2, 
i, j = 1, 2, 3 are defined by 

{7i z/ Re/j(A) < Relj(X) and ReZi(X) > ReZj(X), 
72 if Rek(X) < Relj(X) and Re^(A) < ReZj(X), for A G D n . 
73 z/ Reii(A) > Re^-(A) 

(2.10) 

The following proposition ascertains that the M n 's defined in this 
way have the properties required for the formulation of a Riemann- 
Hilbert problem. 

Proposition 2.1. For each n = 1,2, the function M n (x,t,X) is well- 
defined by equation Ii2.9\) for X G D n and (x,t) G Q. For any fixed 
point (x,t), M n is bounded and analytic as a function of X G D„ away 
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from a possible discrete set of singularities {Xj} at which the Fredholm 
determinant vanishes. Moreover, M n admits a bounded and continuous 
extension to D„ and 



M n (x,t,\)=I + 0(-) 



X — > oo, A G D r , 



(2.H) 



Proof. The bounedness and analyticity properties are established in 
appendix B in [7]. And substituting the expansion 



M = Mn + 



+ 



+ 



A 



oo. 



A ' A 2 

into the Lax pair (12.31) and comparing the terms of the same order of 
A yield the equation (12. 111) . □ 

2.4. The jump matrices. We define spectral functions S n (X), n = 
1, 2, and 

S n (A) = M n (0,0,A), XeD n , n = l,2. (2.12) 

Let M denote the sectionally analytic function on the complex A— plane 
which equals M n for A G D n . Then M n satisfies the jump conditions 



Mi = M 2 J, XeR, 
where the jump matrices J(x, t, A) are defined by 



J 



D (i\Ax+i\Bt) 



(S^Si) 



(2.13) 



(2.14) 



(2.15) 



According to the definition of the 7™, we find that 

/ 73 7i 7i \ / 73 73 7i \ 

7 1 = 73 73 7i 7 2 = 73 73 73 

V 73 73 73 / \ 7i 7i 73 / 

2.5. The adjugated eigenfunctions. We will also need the analytic- 
ity and boundedness properties of the minors of the matrices {/^(x, t, A)}f . 
We recall that the adjugate matrix X A of a 3 x 3 matrix X is defined 

by 

/ mu(X) -m 12 {X) m 13 (X) \ 

-m 21 (X) m 22 (X) -m 23 (X) 

V m 3 i(X) -m 32 (X) m 33 (X) / 



X J 
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where rriij(X) denote the (ij)th minor of X. 

It follows from (J2.3P that the adjugated eigenfunction fj, A satisfies 
the Lax pair 

tf + [i\B 1 pK A ] = -(f r n A . 

where V T denote the transform of a matrix V. Thus, the eigenfunctions 
{nf}\ are solutions of the integral equations 

(4(x,t,X) = I- [ e- ixA ^- x '^ iXB ^' t '\P T dx + Q T )fj, A , j = 1,2,3. 

(2.17) 
Then we can get the following analyticity and boundedness properties: 



fi A : (D 1 ,(H,D 2 ), 

\4 



0, (2.18) 

(I>2,Ml). 



2.6. The computation of jump matrices. Let us define the 3 x 
3— matrix value spectral functions s(X) and S(X) by 

fi 3 (x,t,X) = ^ 2 {x,t,X)e^ XAx+ixM h(X), (2.19a) 

tn(x,t,X) = n 2 (x,t,X)e {iXAx+iX ^S(X). (2.19b) 

Thus, 

s(A)=/i 3 (0,0,A), S(A)=/n(0,0,A). (2.20) 

And we deduce from the properties of fij and \x A that s(X) and S(X) 
have the following boundedness properties: 

s(X): (£>i,0,A), 

S(X): (D^Dt), 

s A (X): (I> 2 ,0,Di), 

S A (A): (Di,0,D 2 ). 

Moreover, 

M n (x,t,X) = ^ 2 (x,t,X)e {ixAx+iX ^S n (X), A e D n . (2.21) 
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Proposition 2.2. The S n can be expressed in terms of the entries of 
s(X) and S(X) as follows: 

( 



So 



m 33 (s)M 2 i (S) -17123 (s) A/31 (S) S13 \ 

„ m 33 (s)Mii(g)-mi3(s)M3i(g) S 23 

621 (s J 'S A )ii (S r s*) 33 

„ m 23 (s)Mii(5)-mi 3 (s)Af 2 i(5) S 33 

i31 (s T S A )n (S r s A )33 

gn mgi (s) A/33 (S)-m 3 i (s) A/ 23 (g) 

(S' J, s^)n (s^Oii 6l3 

S21 m n (s) A/33 (5)-m 3 i (s) A/13 (5) 

0» T S A ) 3 3 



(2.22) 



(S T s A )u 

5*31 



«23 



muMM 23 (S)-m 21 ( S )A/i 3 (S) 



■533 J 



{s T S A ) 33 

where rriij and Mij denote that the (i,j) — th minor of s and S, respec- 
tively. 

Proof. Let 7^° denote the contour (X ,0) — > (x,t) in the (x,t) — plane, 
here X > is a constant. We introduce Hs(x,t, k; X ) as the solution 
of f)2.6p with j = 3 and with the contour 73 replaced by r y 3 °. Similarly, 
we define M n (x, t, A; Xo) as the solution of (12. 9p with 73 replaced by 
7 3 °. We will first derive expression for S n (X; X ) = M n (0, 0, A; X ) in 
terms of S{X) and s(A; X ) = /i 3 (0,0, A; X Q ). Then (T2T22]) will follow 
by taking the limit Xq — > 00. 

First, We have the following relations: 



M n (x,t,X;X ) = fi 1 (x,t. i X)e^ Ax+lX ^R n (X;X ), 
M n {x,t,X;X ) = fi 2 (x,t,X)e^ Ax+a ^S n (X-X ), 
M n {x,t,X;X ) = fi 3 (x,t,X)e^ Ax+tX ^T n (X-X ). 

Then we get R n {X; Xq) and T n (X; Xo) are defined as follows: 



(2.23) 



R n (X; X ) 



iXBT 



M n (0,T,X;X ) 



(2.24a) 



(2.24b) 



T n (A;X ) = e- lXAXo M n (X o ,0,X;X Q ). 
The relations ( I2.23P imply that 

s(X;X ) = S n (X;X )T-\X;X ), S(X) = S n (X; X )R-\X; X ). 



(2.25) 
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These equations constitute a matrix factorization problem which, given 
{s(A), 5(A)} can be solved for the {R n , S n ,T n }. Indeed, the integral 
equations (12.91) together with the definitions of {R n , S n , T n } imply that 

(R n (X;X )) tJ =0 if 7S = 7i, 
(S n (X;X )) tJ =0 if 75 = 72, (2.26) 

(r n (A;X o ))<i = if 7^ = 73. 

It follows that (I2.25P are 18 scalar equations for 18 unknowns. By 
computing the explicit solution of this algebraic system, we find that 
{5 n (A; X Q )}\ are given by the equation obtained from (12.221) by replac- 
ing {5 n (A), s(A)} with {S n (X; X ),s(X;X )}. taking X — > oo in this 
equation, we arrive at (12.221) . □ 

2.7. The residue conditions. Since fi 2 is an entire function, it follows 
from (12.211) that M can only have singularities at the points where the 
S ; n s have singularities. We infer from the explicit formulas (I2.22p that 
the possible singularities of M are as follows: 

• [M]i could have poles in D 2 at the zeros of (5 T s A )n(A); 

• [M] 2 could have poles in D\ at the zeros of (s T 5 A )n(A); 

• [M] 2 could have poles in D 2 at the zeros of (s T S A ) 33 (X); 

• [M] 3 could have poles in D\ at the zeros of (5' T s A )33(A). 

We denote the above possible zeros by {Xj}± and assume they satisfy 
the following assumption. 

Assumption 2.3. We assume that 

• (S T s A )n(X) has uq possible simple zeros in D 2 denoted by 

• (s T S A ) ii (A) has n\ —no possible simple zeros in Di denoted by 

lA'/no+l/ 

• (s T S A )s3(X) has n 2 — ni possible simple zeros in D 2 denoted by 

• (S T s A )33(X) has nj, — n 2 possible simple zeros in D\ denoted by 
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and that none of these zeros coincide. Moreover, we assume that none 
of these functions have zeros on the boundaries of the D„ 's. 

We determine the residue conditions at these zeros in the following: 

Proposition 2.4. Let {M„}\ be the eigenfunctions defined by 112. 9\) and 
assume that the set {Xj}± of singularities are as the above assumption. 
Then the following residue conditions hold: 

A_AjL Ji (5 t s^)ii(a j ) msitv,) L v 3nA ' (2.27a) 
1 <j < n ,\j e D 2 
R es . , \m]o = 1 m 21 (s t s a ) 33 (\j) e 6l2 ^\M(\-)]i 

neS X =X ] [IVl\ 2 -^-^—— (S 13 (A i ) S3 i(A i )-5 33 (A 3 ) S11 (A i )) e L-M^jJJl) 

"o <j < ni,Xj E D 1 

(2.27b) 

Res^[M h = -j^^j^^^ 
n\ < j < n 2 ,\j E D 2 

(2.27c) 

Resx-x[M] 3 = — -J (Sl3a21 - 5 ^ ll)(Aj) e g23 ^)[M(A 7 -)] 2 , 

n 2 < j < N, Xj e D 1 

where f = %., and 9ij is defined by 

6 ij (x,t,X) = (k-l J )x+(z i -z j )t, i,j = 1,2,3. (2.28) 

Proof. We will prove (I2.27al) . (|2.27bl) . the other conditions follow by 
similar arguments. Equation (I2.2ip implies the relation 

Afi = ^ 2 e {aAx+ixiit) S u (2.29a) 

M 2 = fi 2 e ilxAx+lX ^ t) S 2 . (2.29b) 
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In view of the expressions for S\ and S 2 given in (12. 22ft . the three 
columns of (I2.29a|) read: 

[Mi]i = Hi^ii(A) + [^he e21 s 2 i(X) + H3e e31 s 3 i(A), (2.30a) 

[M 1 } 2 = [^ 2 ]ie g " m33 ffi^Ti M31 (A) + l^h ^ps^lf" ( A ) 

;Mii-mi 3 . 
(s T S A )n 

(2.30b) 



_i_r,, 1 r 8:n m 23 M 11 -m 13 M 2 i ( \\ 



[M ll3 = Wl /-^-(A) + [, 2 ] 2e »,^ W + W.(^(A)- 

(2.30c) 
while the three columns of (12.29bl) read: 

[M 2 ]i = Ni77#^(A) + i^he^jj^- (A) + [/i 2 ] 3 e e31 -|V (A), 



(S T s A ) n v v ™ (S¥) u lv "" (S T s 



Ml 



(2.31a) 

[ M 1 _ [,, 1 p 9i2 m 2 i A/33 -m 3 i A/23 f \\ _i_ [.. 1 mi i A/33 -m 3 i A/13 ( \\ 
[M 2 \ 2 - LM2jie (s T 5 A) 33 (AJ + [^ 2 J 2 ( S T 5 A) 33 (AJ 

A/23—?"; 
(s T ^) 3 



1 f„ 1 ^6*32 miiM23-m 2 i A/13 ( \ \ 



(2.31b) 

[M 2 ] 3 = [^}isi 3 e ei3 + N 2 s 23 e e23 + [^W- (2.31c) 

We first suppose that Aj G D2 is a simple zero of (5 ,T s j4 )n(A). Solving 
( I2.31b|) and ( 12.3 left for [// 2 ]i and [/t 2 ] 2 and substituting the result in to 
f l2T3TiD . we find 

Mi = |?gf <^[M 2 ] 2 + Sfe^[M 2 ] 3 

Taking the residue of this equation at Xj, we find the condition (|2.27aj) 
in the case when Xj e D 2 . Similarly, we can get the equation (12.27dj) . 
Then let us consider that Aj G D\ is a simple zero of (s T S A )n(X). 
Solving (I2.30al) and (I2.30cl) for \p2\1 an d [/i 2 ] 3 and substituting the 
result in to (12.30bl) . we find 

M 2 i(S T S A ) 33 D 6i2\MA {S T s A ) 33 



[Mi] 2 = - rr^f s j33 re* u [Mi]i- g ^ % j33 [U2I2 

1 "- (s J 5 A )n(Si 3 S3i-5 3 3Sn) L J-U S13S31-S33S11 L ^ J ^ 

-e*»[Mi] 3 



"t23(-g T S A )33 r 6> 32 [ 
■513531—533511 
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Taking the residue of this equation at Xj, we find the condition ( 12. 2 7b j) 
in the case when Xj e D\. Similarly, we can get the equation (|2.27c[) . 

D 

3. The Riemann-Hilbert problem 

The sectionally analytic function M(x, t, X) defined in section 2 sat- 
isfies a Riemann-Hilbert problem which can be formulated in terms of 
the initial and boundary values of p%j(x, t). By solving this Riemann- 
Hilbert problem, the solution of (11.11) can be recovered for all values of 

x, t. 

Theorem 3.1. Suppose thatpij(x, t) are a solution of ( 11.11) in the half- 
line domain Q with sufficient smoothness and decays as x — > oo. Then 
Pij(x,t) can be reconstructed from the initial value {Pij,o{ x )}i j=i an d 
boundary values {Qij,o(t)}fj =1 defined as follows, 

Pij,o( x ) = Pij( x i 0)> 9ti,o(<) = P»i(0, t). (3.1) 

Use the initial and boundary data to define the jump matrices J(x, t, X) 
as well as the spectral s(X) and S(X) by equation /i2.19\) . Assume 
that the possible zeros {Xj}± of the functions (S' T s A )33(A), (s T S A )n(X), 
(s T S A ) 33 (X), (S T s A ) 33 (X) are as in assumption 2.3. 

Then the solution {pij(x, t)}f J=1 is given by 

Pij(x, t) = —i(di — dj) lim (XM(x, t, X))ij. (3.2) 

A— >oo 

where M(x, t, A) satisfies the following 3x3 matrix Riemann-Hilbert 
problem: 

• M is sectionally meromorphic on the complex X— plane with 
jumps across the contour IR, see Figure 2. 

• Across the contour M., M satisfies the jump condition 

M 1 (x,t,X) = M 2 {x,t,X)J(x,t,X), XeR. (3.3) 

where the jump J is defined by the equation \2.1J$ . 

• M(x,t,X) =I + 0(i), A^oo. 
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The residue condition of M is showed in Proposition\2.4. 



Proof. It only remains to prove ( 13.21) and this equation follows from 
the large A asymptotics of the eigenfunctions. 

We write the large A asymptotics of M as follows: 

M(x, t, A) = M (x, t) + Ml V ' + • • • ■ A^+oo. (3.4) 

A 

And insert this equation into the equation ( 12 .3p and compare the co- 
effients of the same order A, for the 0(A), we have M is a diagonal 
matrix ; for the 0(1), we get Mq = I by comparing the diagonal ele- 
ments, and we can have the following equation by comparing the other 
elements 

-i[A,M 1 \=P, (3.5) 

this equation reads the required result of Pij(x, t) 

Pij(x,t) = -i(a,i - aj)M lyij (x,t). (3.6) 

□ 
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